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I ABSTRACT 

o\ ■ 

. We study the secular evolution and gravitational wave signature of a newly- formed, rapidly ro- 

tating neutron star. The neutron star may arise from core collapse in a massive star or from the 
j-jTj. accretion-induced collapse of a white dwarf. After a brief dynamical phase, the nascent neutron star 

' settles into an axisymmetric, secularly unstable equilibrium state. Gravitational radiation drives the 

star to a nonaxisymmetric, stationary equilibrium configuration via the bar-mode instability. The 
emitted quasi-periodic gravitational waves have a unique signature: the wave frequency sweeps down- 
ward from a few hundred Hertz to zero, while the wave amplitude increases from zero to a maximum 
and then decays back to zero. Such a wave signal could be detected by broad-band gravitational wave 
interferometers currently being constructed. 
^ ' We also characterize two other types of gravitational wave signals that could arise in principle 

from a rapidly rotating, secularly unstable neutron star: a high-frequency (/ ^ 1000 Hz) wave which 
increases the pattern-speed of the star, and a wave that actually increases the angular momentum of 
qq . the star. 

O 
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1. INTRODUCTION 



The advent of the new generation of gravitational wave detectors, such as the Laser Interferometer 
Gravitational Wave Observatory (LIGO; see Abramovici et al. 1992) and its French-Italian counterpart 
VIRGO (Bradaschia et al. 1990), stimulates a renewed interest in astrophysical sources of gravitational 
waves (Thorne 1987, 1992). Of special importance are those sources that radiate in LIGO's sensitivity 
band, 10— 1000 Hz. In this paper we study possible unique gravitational wave signals with frequency 
in this range from newly-formed neutron stars. 

Gravitational collapse leading to the formation of a neutron star has long been considered an 
observable source of gravitational radiation. The event rate of type II supernovae is estimated to be 
several per year out to a distance of 10 Mpc, the distance to the Virgo cluster of galaxies (van den 
Bergh & Tammann 1991). The rate of electromagnetically "quiet" gravitational collapses (such as 
accretion-induced collapse of white dwarfs) may be as high as one per galaxy every one to five years 
(Bahcall & Piran 1983, Blair 1989). However, the strengths of gravitational wave signals from core 
collapse are still highly uncertain because they depend crucially on the unknown asymmetry of the 
collapse. If the collapse is axisymmetric, the energy radiated during the implosion phase of the collapse 
is at most ~ 5x 10~ 8 Mc 2 (Finn 1991, Monchmeyer et al. 1991). Such a supernova would not be 
detectable, even by the advanced LIGO detectors, for sources beyond our galaxy and the Magellanic 
clouds. However, it is widely expected that if the collapsing and bouncing core rotates so rapidly that 
it becomes nonaxisymmetric, the efficiency of producing gravitational wave can be much higher. So 
far, nonaxisymmetric core collapse has not been studied except for the calculations based on classical 
homogeneous ellipsoid models (see Saenz & Shapiro 1981 and references therein). 

It is well known that nonaxisymmetric instabilities can develop in rapidly rotating fluid bodies 
when the ratio /3 = T/|W| of the rotational energy T to the gravitational potential energy W is 
sufficiently large (Chandrasekhar 1969; hereafter referred as Ch69). A star becomes dynamically 
unstable when (3 > (3d yn — 0.27, while its I = m = 2 "bar" mode oscillations become secularly unstable 
due to viscosity (Roberts & Stewartson 1963) or gravitational radiation reaction (Chandrasekhar 
1970a) when /3 > (3 sec ~ 0.14. The gravitational radiation (GR) instability for higher-order modes 
(m > 2) can set in at even smaller value of (i (Friedman & Schutz 1978, Comins 1979a, b). When 
the degenerate core of a massive star or when a white dwarf with mass exceeding the Chandrasekhar 
limit collapses to form a neutron star, the value of /3 can increases significantly. Assuming angular 
momentum conservation, the collapsed star would have 



where Ri and R are the radius of the core before and after collapse, and fli refers to the precollapse 
initial value. Typically, Ri ~ 1500 km for a white dwarf near the Chandrasekhar mass, and R ~ 10 
km for a neutron star; hence Ri/R ~ 10 2 . Thus the nascent neutron star can be rapidly rotating even 
if the initial precollapse core is not. 

While there are no direct observations of the rotation of precollapse cores — they are hidden 
inside the envelopes of giant stars — the situation for the accretion-induced collapse of white dwarf 
and for the collapse of merging white dwarf binaries may be clearer. In both cases, the precollapse 
white dwarfs can attain sufficient angular momentum from the accreting gas or from the initial binary 
orbits. Although most millisecond pulsars are thought to have been spun up by accretion, it has been 
suggested that some millisecond pulsars originate from the collapse of white dwarfs in binary systems 
(Michel 1987, van den Heuvel 1987, Bailyn & Grindlay 1990). The initial spins of ordinary radio 
pulsars are not well constrained by the pulsar statistics; but it is possible that some radio pulsars are 
created with high initial spin. Although the physical mechanism that gives rise to the high velocity 
(<~ 1000 km/s) of some radio pulsars such as PSR 2224+65 (Cordes, Romani & Lundgren 1993) is 
not completely clear (see Radhakrishnan 1992 for a review), it is conceivable that rapid rotation and 
asymmetry of the collapsed core play an important role. Finally, there may already be some indirect 
observational evidence for asymmetric core collapse: the spectropolarimetry of SN1987a (Cropper 
et al. 1988) and SN1993J (Trammell et al. 1993), as well as emission line profiles in the spectra of 
SN1993J (Spyromilio 1994), indicate an asymmetric atmosphere in both supernovae. 
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A new-born neutron star can be secularly unstable but dynamically stable (f3 sec < (3 < Pdyn) 
only if the rotation rate of the precollapse core lies in a narrow range. Thus it is more likely that the 
core becomes dynamically unstable (f3 > (3d yn ) following collapse, provided that the initial is not 
too small. The fate of a dynamically unstable rotating star has been studied numerically by Durisen 
et al (1986), Williams & Tohline (1988), and more recently by Houser, Centrella & Smith (1994), who 
have also calculated the resulting gravitational waves due to the dynamical instability. According to 
Houser et al (1994), the instability develops as follows: after shedding mass and angular momentum 
in the form of two tailing spiral arms, which then expand and merge, causing shock dissipation, the 
system evolves on a dynamical timescale toward a nearly axisymmetric equilibrium state. Its (5 falls 
below the dynamical stability limit f3d yn , but remains much larger than the secular stability limit f3 sec . 
Such a post-collapse, axisymmetric core will evolve into a nonaxisymmetric configuration on a secular 
dissipation timescale. The main purpose of our paper is to study the growth of this secular instability 
(§2) and to characterize the resulting gravitational wave form (§3). We follow the secular evolution 
using a compressible ellipsoid model, which we have recently developed to study the equilibrium and 
dynamics of rotating stars and binary systems (Lai, Rasio & Shapiro 1993, 1994, hereafter referred as 
Paper I and Paper II). We shall see that for the typical neutron star with mass M ~ 1.4M Q and radius 
R ~ 10 km, the gravitational wave amplitude increases from zero to a maximum and then decreases 
to zero again, while the wave frequency sweeps from a few hundred Hertz down to zero (Fig. 4). Such 
a gravitational wave signal can in principle be detected by LIGO (§3.3, especially Fig. 6). 

Another purpose of our paper is to characterize all types of gravitational waveforms one can expect 
from a rotating star modeled as a compressible fluid ellipsoid. The secular evolution of homogeneous, 
incompressible ellipsoidal figures has been studied before (Miller 1974, Detweiler & Lindblom 1977), 
but the waveforms have not been considered or characterized. In addition to the low-frequency wave 
considered in §3, we find two other types of wave forms: a high-frequency wave with decreasing 
amplitude and increasing frequency, and a peculiar wave which increases the angular momentum of 
the star. These are discussed in §4. 

The fluid viscosity can have some effects on the secular instability and the evolution. These effects 
are discussed in §5. A summary including a discussion of the main uncertainties is given in in §6. 

We use units such that G = c = 1 throughout the paper. 

2. GRAVITATIONAL RADIATION INSTABILITY IN 
ROTATING STARS 

Chandrasekhar (1970a) first discovered that GR reaction can induce a nonaxisymmetric secular 
instability in a uniformly rotating, incompressible Maclaurin spheroid. He found that the I = m = 2 
"bar" mode (with angular behavior of the form e ±lm ^) becomes unstable beyond the point at which the 
nonaxisymmetric Jacobi sequence or Dedekind sequence bifurcates from the axisymmetric Maclaurin 
sequence (Ch69). At the bifurcation point, the ratio of the kinetic energy and the absolute value of 
the potential energy (3 = T/\W\ equals f3 sec — 0.1375. 

We have re-analyzed this bar-mode instability induced by GR as well as by viscosity in compress- 
ible rotating stars (Appendix A). Our analysis is based on the equilibrium and dynamical ellipsoid 
model developed in Papers I-II, in which a rotating star is treated as a self-similar ellipsoid with 
polytropic density and linear velocity profiles. We have extended many known results concerning 
incompressible ellipsoidal figures (Ch69) to their compressible analogs. In particular, we find that, 
just as in the incompressible case, the secular instability sets in at the bifurcation point, where 
= Psec = 0.1375 independent of the polytropic index n. Similarly, a dynamical bar-mode instability 
sets in at (3 = (3d yn = 0.2738. While for strictly uniformly rotating stars, mass-shedding can occur be- 
fore the point of bifurcation for n > 0.808 (James 1964; see Tassoul 1978 for review), a slight amount 
of differential rotation can in principle inhibit mass-shedding without changing the global structure 
of the star. Indeed, axisymmetric models of differentially rotating polytropes which retain the global 
properties of Maclaurin spheroids have been constructed by Bodcnheimer & Ostriker (1973). Our 
finding that the secular and dynamical stability limits do not depend sensitively on the polytropic 
index n as well as differential rotation law is also corroborated by earlier numerical studies (Ostriker 
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& Bodcnhcimcr 1973; see also Managan 1985; Imamura, Friedman, & Durisen 1985; Ipscr & Lindblom 
1990). 

Along an equilibrium Maclaurin sequence parametrized by the eccentricity e, the ratio (3 and the 
dimcnsionless angular velocity fiju are given by 
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where g n = (1 — n/5)« n , n n < 1 is a constant ol order unity which depends only on n (see Table 
I in Paper I), and p — 3M/ (AirR 3 ) is the mean density. Equations (2.1)-(2.2) are the same as the 
expressions for incompressible Maclaurin spheroids (Ch69), but they are also valid for compressible 
configurations within our approximation. Because of the centrifugal support, the size of the star 
increases as the rotation rate increases. The mean radius of the rotating star R = (aia2a 3 ) 1 / 3 is given 

by 
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where R a is the radius of the nonrotating star with the same mass. For a given equation of state 
and mass, R is uniquely determined. In a rotating frame with angular frequency f2, our analysis 
(Appendix A) yields the bar-mode frequency 
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and Bu is the index symbol as defined in Ch69 (Chap. 3), and depends only on the axis ratios. 
Gravitational radiation de-stablizes the Dedekind-mode (the mode with plus sign in eq. [2.4]) when 
& 2 M > 2Bn, or, equivalently, when (3 > (3 sec . The growth time tqw for the secular instability is given 

by 
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where M 1A = M/(1.4M ), R w = R o /{10 km), and g = qj 2 G j {ixGp) 1 ' 2 = (4B n - Cl 2 M ) 1/2 
(eq. [2.5]). For < (3 — (3 sec « 1, equation (2.6) can be approximately written as 



TGW-2X 10~ b M^XoCS " PsecY* (s) 



(2.7) 



Figure 1 shows the tqw as a function of /? = T/|W| for n = 0, 0.5, 1. We see that for large 0, 
the growth time tqw is larger for more compressible stars because of the rotation-induced expansion 
of the star (eq. [2.3]). However, for n ~ 0.5 — 1, typical of neutron stars, tqw does not depend 
sensitively on the polytropic index n. Thus for a neutron star with M = 1.4M Q , R a = 10 km, we 
have t gw ~ 7 x 10 4 s for = 0.15, t gw ~ 20 s for /3 = 0.20, and t gw ~ 1 s for (3 = 0.24. 

Higher-order (m > 2) modes can become unstable at a smaller [3. In fact, gravitational radiation 
tends to make all rotating star unstable (Friedman & Schutz 1978). However, the growth time of 
this gravitational radiation instability increases rapidly as m increases (Comins 1979a, b). Moreover, 
viscosity tends to counteract this instability (Lindblom & Detweiler 1977), and the stablizing effect 
becomes more efficient for higher order modes. Thus we will not discuss these high-order instabilities 
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for (3 < [3 sec ; the consequences of these instabilities are not as clear as that of the bar- mode instability. 
Note that for /3 just slightly larger than fi seci the growth time for the m = 2 mode instability is longer 
than that of a m > 2 mode (as seen in Fig. 1, where the result for the m = 3 mode is from Comins 
1979b, for n = 0). Thus our discussion below about the evolution will not be applicable to those very 
small values of (i — f3 sec . However, for /3 > 0.17, the growth time for the bar-mode is much shorter (by 
more than a factor of 10) than that of the m = 3 mode Therefore, for these values of /3, the bar-mode 
instability is much stronger than the higher-order modes, and our discussion of the secular evolution 
is valid. 

Note that we have normalized the mass and radius to those of a cold, nonrotating neutron 
star. Following its implosion, which takes ~ 1 s, the 1.4M© collapsing core settles into hydrostatic 
equilibrium within a few milliseconds. This hot protoneutron star does not resemble a usual cold 
neutron star, as it consists of an unshocked cold inner core with radius <~ 10 km surrounded by a 
shocked hot outer core (mantle) which extends to ~ 100 km. The protoneutron star evolves quasi- 
statically via neutrino cooling. However, within one second, the outer mantle cools and contracts 
significantly, and the protoneutron star (with mass ~ 1.4M Q ) approaches a canonical radius of 10 
km. By this time, the hydrostatic structure of the neutron star more or less resembles that of a cold 
neutron star, although the star will remain very hot (T <; 10 10 K) for more than a few tens of seconds 
(Burrows & Lattimer 1986, Burrows 1990). Since the nonaxisymmetric instability growth time is 
typically larger than a second, it is appropriate to use the canonical radius of a cold neutron star (but 
see §6). 

3. EVOLUTION AND GRAVITATIONAL WAVEFORM 

Here we discuss the evolution of a rotating star modeled as a non-axisymmetric ellipsoid (i.e., a 
Riemann-S ellipsoid in the language of Ch69 and Paper I) under the influence of gravitational radiation 
reaction and zero viscosity. 

3.1 Secular Evolution 

A secularly unstable Maclaurin spheroid will evolve away from the axisymmetric configuration 
due to gravitational radiation. The evolution can in principle be studied using the full dynamical 
equations of ellipsoidal figures (Ch69, Miller 1974, Paper II) including gravitational radiation reaction. 
However, since the growth time of the instability tqw is generally much longer than the dynamical 
time of the star, the evolution is quasi-static, i.e., the star will evolve along an equilibrium sequence of 
Riemann-S ellipsoids (Ch69, Paper I). A general Ricmann-S ellipsoid is characterized by the angular 
velocity Jl = fle z of the ellipsoidal figure (the pattern speed) about a principal axis and the internal 
motion of the fluid with uniform vorticity £ — £e 2 along the same axis (in the frame corotating with 
the figure). The vorticity in the inertial frame is then = 2Q + (. Since the gravitational radiation 
reaction acts like a potential force, the fluid circulation C along the equator of the star is conserved 
in the absence of viscosity (Miller 1974; Paper II). Following Paper I (§5.1) we write 

C = ^--L Kn M^ C = (--^-KnM^ 7ra ia2 C (0) =IA- ^K n M ai a 2 n, (3.1) 

where / = n n M{a\ + a\)/5 is the moment of inertia, and A = — aia 2 C/(«i + The quantity C has 
the dimension of angular momentum but is proportional to the conserved circulation C = ira\a2C^ ■ 
We refer to C itself as the circulation for convenience. Note that in the same notation, the angular 
momentum of the ellipsoid is 

2 

J = IQ,— -K n M ai a 2 A. (3.2) 



Thus, under the influence of gravitational radiation reaction, and when the viscous dissipation is 
negligible, a secularly unstable Maclaurin spheroid will evolve along a constant-C equilibrium sequence, 
and proceed ultimately toward a Dedekind ellipsoid, for which Q = (a "stationary football"). The 
evolution track can be best illustrated by the diagram shown in Figure 2, which depicts the equilibrium 
energies for various constant— C sequences as functions of the axis ratio c^/ai, all for n = 1. For 
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comparison, the Jacobi and the Dedekind sequences are also shown, while the Maclaurin sequence 
corresponds to the vertical line with a^/ai = 1. A given Maclaurin spheroid has a unique value of C, 
but two different constant— C sequences branch off: one sequence is Jacobi-like, with \(\ < 2|f2|, the 
other is Dedekind-like, with |£| > 2|f2|. Clearly, the endpoint of the evolution depends on the value \C\ 
of the initial configuration. For \C\ < \C sec \, where C sec is the circulation of a Maclaurin spheroid at 
the bifurcation point, the Maclaurin spheroid is the lowest energy terminal state toward which both 
Jacobi- like and Dedekind like configurations evolve; for \C\ > \C sec \, on the other hand, the terminal 
state is a Dedekind ellipsoid. A Jacobi-like configuration with \C\ > \C sec \ first evolves toward an 
axisymmetric, Maclaurin state, and then continues through a Dedekind-like phase toward a Dedekind 
ellipsoid. The value of |C sec | = J sec is given by 

J sec = 0.304 ( -^-) V2 (M 3 ^) 1 / 2 , (3.3) 



l-n/5 

(see eq. [3.27] and Table II in Paper I), where R is given by equation (2.3). 

Now let us focus on the evolution of a secularly unstable Maclaurin spheroid toward a Dedekind 
ellipsoid, corresponding to region I in Figure 2. The initial Maclaurin spheroid has angular momentum 
Ji = \&\. The final Dedekind state has J f = \C t \[2 ai a 2 / \a\ + a 2 .)]/ (cf. cqs. [3.1]-[3.2] with Q, = 0), 
the total angular momentum radiated in gravitational wave is given by 



A J = Ji 



(ai - a 2 ) 2 



(3.4) 

/ 



where the subscript "f" stands for the final Dedekind configuration. The total energy AE emitted 
as gravitational radiational wave for a secularly unstable Maclaurin spheroid can be calculated by 
comparing the difference in equilibrium energy between the Maclaurin spheroid and the Dedekind 
ellipsoid with the same C. Figure 3 shows the dependence of AE as a function of [3 = T/\W\ of the 
initial Maclaurin spheroid. We see that AE does not depend sensitively the compressibility for n <J 1. 
Thus the energy radiated in gravitational wave during the Maclaurin-Dedekind transition can be as 
large as 0.02M 2 /R o , or AE <, 4 x 10~ 3 M for a typical neutron star radius (R /M ~ 5). This is much 
larger than the energy radiated (<, 5 x 10~ 8 M) in a axisymmetric collapse preceding the neutron star 
formation (Finn 1991, Monchmeyer et al. 1991). 

3.2 Gravitational Wave Amplitude 

The gravitational wave emitted during the Maclaurin-Dedekind evolution has certain unique 
characteristics. At the beginning of the evolution, the wave amplitude is small since the star is still 
nearly axisymmetric. Near the end of the evolution, as the star approaches a Dedekind configuration 
asymptotically, the wave amplitude is small again because the angular velocity f2 of the figure becomes 
decreasingly small. Thus we expect a non-monotonic behavior of the gravitational wave amplitude 
during such a evolution, while the wave frequency decreases monotonically. 

In the weak-field, slow- motion limit, the gravitational radiation can be calculated by the standard 
quadrupole formula (see, e.g., Misner, Thorne & Wheeler 1970). The rates of energy loss and angular 
momentum loss are given by 

(SL-(SL 

where la = n n Ma 2 /5 are the components of the star's quadrupole moment along the principal axes 
in its equatorial plane. At a distance D from the source, the two polarization components of the 
waveforms are 

h+ - — |fi 2 (/n -I 22 )cos$(l + cos 2 60, 

4 (3-6) 
h x = -— f2 2 (/n - I22) sin$cos#, 
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where is the angle between the rotation axis of the star and line of sight from the earth, and 
$ = 2 J VLdt is twice the orbital phase. The wave is quasi-periodic, with frequency / = Q/ir and 
amplitude (for 8 = 0) 

h=^n 2 (I 11 -h 2 ). (3.7) 

Typical wave amplitudes as a function of frequency are shown in Figure 4, for polytropic index 
n = 0.5 and 1. The initial Maclaurin spheroids have (3 = 0.2 and (3 = 0.24, respectively. The wave 
sweeps with time from high frequency to low, that is, from right to left in Figure 4. Also shown in the 
figure are the timescale of the evolution as measured by \dt/d\nf\ = |///|, and the number of cycles 
gravitational wave spent near frequency / 
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(3.8) 



where dE/dfl is evaluated along an equilibrium sequence with constant C. The number of cycles near 
the frequency at which the wave amplitude is maximum depends sensitively on the initial value of [3, 
and to a less extent, on the polytropic index. But even for the initial (3 close to the dynamical limit 
(3dyni there are still 10 3 cycles when the wave amplitude is close to the maximum. Note that at the 
beginning of the evolution (/ — > f m ax), the timescale approaches infinity. This is because we have 
assumed exact equilibrium models, and a perfect Maclaurin spheroid does not evolve at all. In reality, 
a small perturbation will make the Maclaurin spheroid evolve on the order of a few growth time tqw , 
as determined in §2. This time is of the same order of magnitude as the time the system spends near 
the maximum amplitude phase. 

It is important to note that the wave frequency is determined by the pattern speed Q of the 
ellipsoidal figure, and this frequency can be much lower than the rotational frequency of the secularly 
unstable star. During the evolution, the wave frequency decreases monotonically from an initial 
value to zero. Consider the maximum frequency of the wave, corresponding to the frequency near 
the beginning of the evolution, when the star is still close to a Maclaurin spheroid. The bar-mode 
oscillation frequency as viewed in a general rotating frame is given by equation (2.4), and gravitational 
radiation destablizes the Dedckind-mode. In a particular reference frame, with angular velocity £l ma x, 
the Dedckind-mode can be neutralized, i.e., a+(fl ma x) = 0. It is along the axes of this rotating 
frame that the Maclaurin spheroid starts deforming into a nonaxisymmetric shape. Therefore, as the 
Maclaurin spheroid begins to evolve due to gravitational radiation, the pattern speed Q m ax is given 

by 



£7 r t 



-{Q.M — Co), 



(3.9) 



which is always smaller than the initial rotation rate Qm of the star. Clearly, at the bifurcation point, 
Qmax = (i.e., at the bifurcation point, the oscillation mode is neutralized in the inertial frame). The 
angular frequency Q ma x is also one of the two values assigned to a Maclaurin spheroid if it is to be 
considered as a limiting member of Riemann ellipsoids (Ch69, §48). Equation (3.9) is simply half of 
the absolute value of the mode frequency a + (0) in the inertial frame. The maximum frequency of the 
gravitational wave during the Maclaurin-Dedekind evolution is then 
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Figure 5 shows f max as a function of initial (3 for n = 0, 0.5 and 1. For a given [3, the maximum 
frequency is smaller for larger n, as a result of the rotational expansion factor R/R Q . This dependence 
on n is significant only for larger f3 and n <; 1. For comparison, we also show in Figure 5 twice of 
the rotational frequency of the Maclaurin spheroid 2/m = ^m/tt- This would be the frequency of the 
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gravitational wave emitted if the rapidly rotating star were deformed slightly away from axisymmetry 
by external stress (e.g., a mountain on the neutron star). This frequency is larger than 1000 Hz if 
star rotates sufficiently rapidly to be secularly unstable. 

3.3 Characteristic Wave Amplitude 

Thorne (1987) has emphasized the importance of the signal-to-noise ratio of a given source when 
considering its detectability. For a broad-band detector such as LIGO, the best signal-to- noise ratio 
will be obtained by matched filtering of the data, S/N ~ h c /h rmsi where h rms is the interferometer's 
root mean square noise (Abramovici et al 1992), and h c is the characteristic amplitude of the source. 
For the quasi-static evolution discussed in §3.2, h c is given by 3 

1/2 - - (— — y /2 - — /^v / Y— — v /2 

C ~ d\nf ~ D\2k dtl eq ) ~ D\M) \2n dti e J ( ' ' 

where E = E/(M 2 /R ), fl = ^/(M/i? 3 ) 1 / 2 , and we have used equations (3.7)-(3.8). 

In Figure 6, we show h c for the four wave amplitudes depicted in Figure 4, and compare them 
with the rms noise h rms and the sensitivity to bursts h s b — llh rms of the advanced LIGO detec- 
tor (Abramovici et al. 1992). From our numerical calculations, we find (dE/dri) eq oc £l to a good 
approximation. Thus the characteristic amplitude of gravitational wave during the evolution from a 
Maclaurin spheroid to a Dedekind Ellipsoid is given by 

h c * 6.0 x 10- 3 (22^2) Mf 4 4 < 4 / 1/2 , (3.12) 

where / is the wave frequency in Hertz. This expression is accurate to within ~ 20% to all relevant 
values of (3 and n (0.16 < (3 < 0.27, < n < 1.5). We see clearly from Figure 6 that the gravitational 
waves from the non-axisymmetric evolution of rapidly rotating neutron stars should be detectable with 
high confidence (S/N > 11) out to the distance of ~ 140 Mpc by the advanced LIGO detector at ~ 100 
Hz. Valuable information about the masses and radii, as well as the rotation rates of newly-formed 
neutron stars may be gained from the unique wave signatures (Figs. 4-5). 

4. OTHER TYPES OF WAVEFORM 

The low-frequency gravitational wave discussed in §3 results from the Maclaurin-Dedekind transi- 
tion (region I in Fig. 2). Apart from this low- frequency wave, there are two other types of gravitational 
waveforms that one may expect, at least in principle, during the evolution of a Riemann-S ellipsoid 
due to gravitational radiation. 

• "Spin-up" Waves: A Jacobi-likc ellipsoid (|£| < 2|f2|) will evolve to become a stable Maclaurin 
spheroid if \C\ < \C sec \. If \C\ > |C sec |, it will evolve toward a secularly unstable Maclaurin spheroid, 
and then, after a nonaxisymmctric perturbation, it will further evolve to become a Dedekind ellipsoid. 
The evolution of a Jacobi-like ellipsoid toward a Maclaurin spheroid (regions II and III in Fig. 2) is 
accompanied by emitting gravitational wave with increasingly high frequency. 

A special case is the evolution of a Jacobi ellipsoid. From Figure 2, we see that a Jacobi ellipsoid 
will evolve toward a secularly stable Maclaurin spheroid 4 . Chandrasckhar (1970b) first noted that 
the rotational frequency increases as the Jacobi ellipsoid evolves under gravitational radiation. This 
spin-up is caused by the decrease of the momentum of inertia as the angular momentum is radiated 
away. However, it was assumed that the evolution is along the equilibrium Jacobi sequence. Thus 

3 More precisely, one needs to average over the source direction to obtain h c (Thorne 1987), but 
this only introduces a factor ^/4/5 in eq. (3.11). 

4 For a more compressible configuration with n ^ 1.5, a Jacobi ellipsoid will first evolve toward a 
secularly unstable Maclaurin. This dependence on compressibility comes from the expansion of the 
star as it gets more deformed due to rotation or internal motion. See Fig. 2 in Paper II. 
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the evolution ends up at the bifurcation point. This assumption is not valid, since C is not conserved 
along a Jacobi sequence. The correct evolutionary sequence was computed by Miller (1974) for 
incompressible configurations. 

Ipser & Managan (1984) also considered gravitational wave emission from the evolution of a 
compressible Jacobi ellipsoid, based on their numerical models of rotating polytropes. They too 
made the assumption that the evolution is along a Jacobi sequence, terminating at the bifurcation 
point. Because of the mass-shedding constraint, exactly uniform rotating Jacobi ellipsoid exist only 
for n <; 0.8. Even for n = 0.5, the initial uniformly rotating Jacobi ellipsoid is very close to the 
bifurcation point. Thus Ipser & Managan (1984) found that only a small fraction (~ 10~ 4 — 10~ 3 ) 
of the mass energy can be radiated as gravitational wave, and the wave is nearly monochromatic 
with / ~ fl sec /n <~ 2300 — 2500 Hz (see Fig. 5(b)). This may be somewhat misleading. Indeed, if 
slight differential rotation can inhibit the mass shedding of a Jacobi ellipsoid, as argued in §2, a more 
deformed triaxial Jacobi configuration might exist in principle. As such a highly deformed Jacobi 
ellipsoid evolves toward a Maclaurin spheroid, much more energy can be emitted, and the frequency 
can span a wider range. 

In Figure 7 we show the gravitational waveforms when a Jacobi-like triaxial ellipsoid evolves 
toward a Maclaurin spheroid for n — 1. The terminal Maclaurin spheroids (with the same C as the 
initial states) have (5 = 0.02, 0.12, 0.135 and 0.24. In the first three cases, the evolution ceases once 
the Maclaurin branch is reached , while in the last case, the Maclaurin spheroid (with (3 = 0.24) can 
evolve further to a Dedekind ellipsoid (see §3). In the cases with (3 = 0.12 and 0.135, the initial state 
is chosen to be a Jacobi ellipsoid. The evolution is from low-frequency to high (from left to right in 
Fig. 7). This is very different from the cases discussed in §3.2 (see Fig. 4). The terminal frequency 
ft is determined by the pattern speed fl t of the ellipsoid as it approaches a Maclaurin shape. This 
pattern speed is equal to the angular frequency of the rotating frame in which the Jacobi mode (with 
— sign in cq. [2.4]) becomes neutralized, i.e., a -(fit) = 0. Thus we have 



and the terminal (maximum) wave frequency is ft = fit /ft- Figure 8 depicts ft as a function of of 
the Maclaurin spheroid. Note that at the bifurcation point, we have fit = flui = flsea an d at the 
dynamical stability limit (<r = 0), we obtain O t = fl max = Qm/2. 

As seen from Figure 7 and Figure 4, the evolution of a Jacobi-like ellipsoid toward a Maclaurin 
spheroid is generally much faster than the evolution of an secularly unstable Maclaurin toward a 
Dedekind ellipsoid. This is because a typical "Jacobi-like" configuration has larger \ fl\ than a typical 
"Dedekind-like" configuration, and because the energy radiation rate (eq. [3.5]) depends on fl through 
a higher-power law. For the same reasons, the typical wave amplitude during the "Jacobi-like" evo- 
lution is typically more than ten times larger than the wave amplitude during the "Dedekind-like" 
evolution. However, the frequencies of the wave are much larger for the "Jacobi-like" evolution, 
typically / £ 1000 Hz. 

The characteristic wave amplitude h c during the "Jacobi-like" evolution discussed here can also 
be calculated using equation (3.11). These are also plotted in Figure 6. Again, h c can be fitted to the 
form 



This expression is approximately valid for all values of C and n, and the differences resulting from 
different C, n and initial conditions lie in the different wave frequency band. Although h c during a 
"Jacobi-like" evolution is larger than h c during a "Dedekind-like" phase, the later is actually easier 
to detect since frequency lies around 100 Hz, at which LIGO detectors are most sensitive. 

• Waves that Carry Negative Angular Momentum: The third type of gravitational wave results 
from the evolution in region IV of Figure 2. The initial configuration is Dedekind-like (|£| > 2|f2|), 
but with the pattern speed fl < 0, i.e., fl is opposite to the internal vorticity £ which provides 
major contribution to the the angular momentum J of the system. Such a configuration will evolve 
toward a Dedekind ellipsoid (with terminal fl = 0) when \C\ > \C sec \, or toward a stable Maclaurin 



(4.1) 




(4.2) 
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spheroid (with terminal Q < 0; see eq. [3.9]) when \C\ < \C sec \. Figure 9 depicts two examples of such 
evolution. Note that does not change monotonically as the 02/01 decreases. For an extremely large 
deformation (small 02/01), J is in the same direction as f2, and \ J\ decreases and \Q\ increases during 
the evolution. When the deformation is not too large, the total angular momentum of the system 
is positive due to the large positive internal vorticity (. As the system evolves toward the terminal 
state (Dedekind or Maclaurin), the energy decreases but the angular momentum \J\ increases, while 
the ellipsoid spins down (|0| decreases). This peculiar behavior comes about because fl is negative, 
and dE = SldJ must be satisfied for secular evolution via gravitational radiation dissipation (Ostriker 
& Gunn 1969). The gravitational wave carries away negative angular momentum, resulting in an 
increase of the angular momentum of the system. The wave frequency can also be low, but the wave 
amplitude decreases monotonically. 

It is not clear whether the waveforms discussed in this section is ever relevant to a new-born 
neutron star. This depends on whether the collapsed core is axisymmetric or nonaxisymmctric after 
hydrostatic equilibrium is established following the initial implosion and/or the dynamical instability. 
Earlier calculations by Duriscn et al. (1986), Williams & Tohlinc (1988) found that a dynamically 
unstable rotating star evolves toward a triaxial bar surrounded by a ring. Recent calculations by 
Houser et al (1994), which include shock dissipation of the ejected gas, suggest that the final outcome 
is an axisymmetric core surround by a shocked halo. Clearly this issue will be resolved only by three- 
dimensional collapsing simulations, which have yet to be performed 5 . Nevertheless, the gravitational 
waveforms discussed in §3.3 and in the present section represent the three classes of waveforms that 
can be expected from the secular non-axisymmetric evolution of any rotating star. 



So far in our discussion we have assumed that the star consists of inviscid fluid. When the fluid 
viscosity is sufficiently large, it tends to counteract the GR instability (Lindblom & Detweiler 1977). 
From analysis in Appendix A, the combined effect of viscosity and GR reaction is determined by 
the ratio of the gravitational radiation timescale and the viscous timescale, Q ~ tQ\y/t visc , or, more 
precisely (eq. [A33]) 



where v is the mean kinematic shear viscosity in units of cm 2 /s. Figure 10 shows the critical T/lTF] 
for instability as a function of Q for a Maclaurin spheroid with n = 0, 0.5, 1 and 1.5. For Q <, 1, the 
instability is mainly caused by GR reaction in the "Dedekind- like" mode, while for Q <; 1 , it is mainly 
induced by viscous stress in the "Jacobi-likc" mode. In the limits of Q >> 1 and Q « 1, both 
instabilities set in at the bifurcation point f3 sec = 0.1375. For Q ~ 1, both modes are stablized, and 
the Maclaurin spheroid can be stable all the way to the dynamical limit for certain value of Q <~ 1. 

The secular evolution discussed in §2-4 typically occurs during the first 1 — 100 seconds after 
the neutron star's birth. During this epoch, the neutron star remains very hot, with temperature 
T 10 10 K (Burrows & Lattimer 1986). The shear viscosity resulting from neutron-neutron scattering 
is v n ~ 10p^ 4 T 10 2 cm 2 /s (Flowers & Itoh 1976, with the fitting formula from Cutler & Lindblom 1987), 
where P15 is the density in units of 10 15 g/cm 3 , and T10 = T/(10 10 K). This is much too small to 
have any effect on the secular instability and evolution discussed in §2-4. The neutrino-induced shear 
viscosity (resulting from neutrino-neutron scattering) v v ~ 10 6 pi 5 T10 cm 2 /s (Thompson & Duncan 
1993) is also too small to be of any importance. 

There have been some discussions about the importance of neutron star bulk viscosity in the 
literature. For a hot neutron star, a bulk viscosity can arise from the phase lag between density and 
pressure perturbation due to the relatively long time scale of weak interactions to re-establish chemical 
equilibrium (Saywer 1989). In our ellipsoid model, the bar-mode is associated with the perturbations 

5 In the binary coalescence calculation of Rasio & Shapiro (1994), a Jacobi-like triaxial system is 
produced for a sufficiently small polytropic index. However, this calculation starts from a corotating 
binary near contact. Calculations starting from non-synchronized spins (more realistic for neutron 
stars) may yield a different result; see Shibata, Nakamura & Oohara 1992. 
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of the three axes given by Sa\ = —80,2, 60,3 = 0. Therefore the bulk viscosity has no effect on the 
damping/growth time of the mode. In realistic calculations, however, there is a small but nonzero 
compression/expansion of the fluid associated with the bar mode, thus bulk viscosity can play a role 
(fpscr & Lindblom 1991). Using the Sawyer's formula for bulk viscosity vb ~ 10 1 Vis^io cm 2 /s (where 
we have taken the perturbation timescale to be of order a millisecond), Ipser & Lindblom (1991) found 
that the bulk viscosity can suppress the GR driven instability for T <; 2 x 10 10 K 6 . However, it should 
be noted that the above formula for bulk viscosity assumes that the neutrinos produced in electron- 
capture escape the neutron star freely. Thus the formula is valid only for relatively low temperatures 
(less than a few times 10 9 K; e.g., Shapiro & Teukolsky 1983). During the early epoch of the neutron 
star studied in this paper, the neutrino optical depth is large, and the bulk viscosity is likely to be 
greatly suppressed due the blocking of the neutrino phase space. A more detailed study is needed to 
fully address this issue. 

In the presence of a shear viscosity, the fluid circulation evolves according to 



(cf. Paper II). This result and equation (3.5) for dJ/dt govern the secular evolution of a general 
ellipsoid under the combined effects of viscosity and GR reaction. Some evolutionary tracks have 
been considered by Detweilcr & Lindblom (1977) for incompressible configurations. 

Now let us consider the effect of a small viscosity on the secular evolution. Without viscosity, we 
see from §3 that a Dedckind ellipsoid is one of the possible final states of the evolution. When viscosity 
is also present, however small, only a secularly stable Maclaurin spheroid can be the final state — it is 
the only configuration which does not radiate GR or dissipate energy viscously. As the star approaches 
a Dedekind ellipsoid, the gravitational evolution timescale increases. When this timescale becomes 
comparable to the viscous dissipation timescale, the star will be driven along a nearly-Dedekind 
sequence, eventually toward the bifurcation point. Thus we can expect that many rapidly rotating 
nascent neutron stars eventually settle down at the bifurcation point along the Maclaurin sequence. 
Of course, the stars can then be subjected to further spin-down due to magnetic braking. 



Thorne (1987, 1992) has summarized the three radically different scenarios of rapidly rotating 
collapse and the resulting wave characteristics: (i) the star remains axisymmetric throughout the 
collapse, resulting in a weak wave signal; (ii) the star becomes a (Jacobi-like) triaxial object, emitting 
nearly monochromatic waves with / <~ 1000 Hz; (iii) the core breaks up into lumps, resulting in the 
strongest wave signal. Regarding (ii), we have shown in §4 that in general the wave emitted by a 
Jacobi-like rotating star can span a wide frequency band up to a few thousands Hertz. Moreover, 
we have considered in this paper the gravitational wave signal in another likely scenario for the post- 
collapse, quasi-static stage: After the hydrodynamical collapse which lasts ~ 1 second, the core settles 
down into a nearly axisymmetric, equilibrium state, which is secularly unstable to nonaxisymmetric 
perturbations. Prior to this, the star may have gone through a nonaxisymmetric dynamical phase, 
but the development of the dynamical instability can only render the final configuration dynamically 
stable and not necessarily secularly stable (Houser et al 1994). Due to the bar- mode instability of 
gravitational radiation reaction (Chandrasekhar 1970a, Friedman & Schutz 1978), the nascent neutron 
star evolves quasi-statically toward a Dedekind-like object on a timescale of tens of seconds to a few 
minutes, emitting gravitational waves with frequency sweeping from a few hundred Hertz down to zero, 
and a non-monotonic wave amplitude. The viscous forces finally bring the Dedekind-like neutron star 
to an axisymmetric Maclaurin spheroid near the bifurcation point on a much longer, viscous timescale. 

Gravitational radiation resulting from the Chandrasekhar-Friedman-Schutz instability has been 
discussed in the context of accreting neutron stars (Wagoner 1984). As a slowly rotating neutron star 

6 Ipser & Lindblom only considered uniformly rotating models and were thus restricted to rotation 
frequencies below the mass shedding limit. 
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accretes mass from a binary companion, it can be spun up until it becomes rotationally unstable on 
an accretion timescale. A monochromatic gravitational wave is then radiated, carrying off the angular 
momentum of the accreted material. During accretion, which increases the star's angular momentum 
very slowly, the neutron star can evolve only slightly above the stability limit. Accordingly it is often 
the instability of a higher-order (m = 4 — 6) mode (which sets in at smaller f3) that is relevant, and 
the star can never reach the bar-mode (m = 2) stability limit (at (3 = [3 sec ). Much effort has been 
devoted in calculating accurately the stability limits of these high-order modes (e.g., Ipser & Lindblom 
1990, 1991). By contrast, in the core collapse considered in this paper, the nascent neutron star can 
easily reach a value of /3 much larger than the secular stability limit (3 sec , as a result of the much 
violent initial implosion. Thus in this case, the dominant instability is the bar-mode. Schutz (1989) 
has also noted that the instability-driven spin-down can yield low frequency gravitational waves, but 
he has focused on the high-order modes. As we noted before, unlike the bar-mode instability, the 
quantitative outcome of a high-order instability is not yet clear at present. 

One important question we have not addressed so far is the size R of the collapsed core. The 
gravitational wave frequency depends sensitively on this quantity through / oc Rq 3 ^ 2 (cf. cq. [3.10]). 
We have adopted R a = 10 km throughout the paper. While the thermal pressure is unlikely to affect 
the radius of the collapsed core appreciably (§2), it is possible that the collapse can be stopped by 
centrifugal forces at sub-nuclear density. Such "centrifugal" bounce has been found in some numerical 
simulations of axisymmetric collapse (Monchmcyer et al. 1991, Yamada & Sato 1993). If this indeed 
occurs, the wave frequency would be lower and the evolution timescale longer than discussed in this 
paper, although the characteristics of the waves do not change. On the other hand, it is conceivable 
that any realistic rotating core will ultimately contract to typical neutron star size. This must be 
the case if we believe that most neutron stars are born with high spin rate. Actual detection of 
gravitational waves from newly-formed neutron stars might provide valuable information about the 
size and rotational properties of the collapsed cores and the dynamics of core collapse itself. 
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APPENDIX A: THE SECULAR INSTABILITY GROWTH TIMES OF 
COMPRESSIBLE MACLAURIN SPHEROIDS 

To calculate the oscillation modes of a Maclaurin spheroid, we start from the general dynamical 
equations for a Riemann-S ellipsoid, including the viscous forces and gravitational radiation (GR) 
reaction (see Paper II, §2,4): 7 
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7 A bulk viscosity contribution, which was not included in Paper II, can be easily incorporated by 
adding a term — (5Db/ K n )(ai/ai + d 2 /a2 + dz/az)/a\ to the right hand side of eq. (Al) and similar 
terms to eqs. (A2)-(A3), where vb is the mean kinematic bulk viscosity. However, this will have no 
effect on the bar-mode in which we are most interested here. 
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The dynamical varibles are the three axes a, (i = 1,2,3), the angular frequency of the figure (pattern 
speed) 0, and the angular frequency A specifying the internal fluid motion (relative to the figure) 
with uniform vorticity £ = — \{a\ + a 2 )/aia 2 ]A. The other quantities have the same meaning as in 
Paper II. 

When linearizing the dynamical equations, we must allow for the fact that a Maclaurin spheroid 
need not be assigned the angular frequency of the frame in which it is at rest. When viewed from 
a rotating frame with angular frequency f2 relative to the inerial frame, the Maclaurin spheroid 
will appear as having stationary internal motion with the uniform vorticity ( — 2(0 jvj — ty, i- e -> 
A = Q — fi M , where U M is the angular frequency of the Maclaurin spheroid relative to an inertial 
frame. Thus the values of Vt eq and A eq for the equilibrium Maclaurin spheroid can in principle be 
arbitrary, subject only to the condition 



^eq A eq — Ojy. 



Consider small perturbations of the form 



a-i = a ieg (l + on), tt = tt eq +uj, A = A eq + \ 
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with en oc e~ lat and similarly for u>, A. Note that tt eq (or A eq ) must be determined from the linear 
mode analysis. We are only interested in the non-axisymmetric bar-mode (called "toroidal mode" in 
Ch69) for which it turns out a\ — — a 2 and a 3 = (see Paper II, §3.3). 

A.l Pure Dynamical Oscillations 

First consider the case when viscosity and GR reaction are absent. A useful expression is 

Oj dAi 
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(no summation for repeated indices), where Ai, Bij are defined as in Ch69 (Chap. 3). Also it is useful to 
note that the pressure term in (A1)-(A3) can be written as hk\P c / '(nn n p c ) — (M/q n R )(R / R) 3 / n — 
5U/(nK n M), where U is the internal energy (Paper II, eq. [2.9]). Linearizing equation (AI) and using 
the equilibrium conditions (Paper II, cq. [3.24]), we obtain 
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where the notation is the same as in Paper II (§3.3), with In — n n Ma\/5 and £ = — 3M 2 /[2(5 — n)B?\. 
The equation for a 2 can be obtained by inter-changing subscript 1 and 2 in (A9). Linearizing equations 
(A4)-(A5) yields 

V eq = -A eq = ^f, (A10) 
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and 

u- A = -n M (ai +a 2 ), 

where we have used f2 eg — A e? = Qm (eq. [A6]). Substituting (AlO)-(All) into (A9) and the corre- 
sponding equation for a 2 , we can see that ct\ = —a 2 , a 3 = is indeed a solution. After some algebra, 
we obtain 

JnAS - -hiol&a, {A12) 



where Aa = a\ — a 2 , and 
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Thus the mode frequency is a — ±ct , in agreement with the result of Paper II, obtained using a 
Hamiltonian formalism. Note that the frequency a is the viewed in a particular rotating frame with 
il = fijvf/2, in which the oscillation is self-adjointed (Cl — —A), and the directions of the principal 
axes of the ellipsoid are fixed (Ch69, §36). In a frame with a different 0, the oscillation frequency is 



since this mode corresponds to a I = m = 2 mode with a exp(z20) dependence on the azimuthal angle 
(see Ch69, §36; Rossner 1967). 

A. 2 Effect of Viscosity 

Now consider a small viscosity v so that we only consider the leading correction ~ 0(p) to 
the mode frequency. Again assuming a e~ l<yt time dependence of the perturbations, from equations 
(A4)-(A6) we obtain 1 

Qeg = (l + i 5 \ ) , A eg = -^0 M (l-i ^2 ), (-415) 

while equation (All) still holds. The fact that Q eq contains an imaginary part implies that the 
rotation of the frame in which the oscillation is specified can itself enhance or damp the oscillation 
amplitude. Substituting (All) and (A15) into the linearized equations for d\, we obtain 

Ad = -o 2 Aa - -^Ad, (A16) 

from which we have 

5u 

a = ±a -i ( A17 ) 

However, this is the frequency of oscillation in a frame with angular frequency given by a complex Q, eq . 
Similar to equation (A14), in a general frame with (real) angular velocity f2, the oscillation frequency 
is 
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Clearly, the Dedekind-mode (with upper sign in eq. [A18]) is always secularly stable. For the Jacobi- 
mode (with lower sign in eq. [A18]), secular instability sets in when Qm > ct , or when 

q j3l > 2B 1U (A19) 
up 
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corresponding to bifurcation values e = 0.8127 and T/|W| = 0.1375. The instability growth time t v 
is given by 



hv ( VLm — o 
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A. 3 Effect of Gravitational Radiation 

(5) 

Now to include the GR reaction, we first need to calculate f ^ , the fifth time derivative of the 
reduced quadrupole moment tensor of the body in the inertial frame projected onto the rotating body 
frame (with basis vectors along the principal axes). This is done using the procedure given in Miller 
(1974) and Appendix A of Paper II. Unlike the quasi-static case (as in Appendix A of Paper II), in 
order to discern the instability, 1 \- need to be calculated to all order. For small perturbations around 
a Maclaurin spheroid, we obtain 
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where on the right-hand side, the superscript (i) now means i-th time derivative in the frame of the 
body. Thus the only non-trivial components are therefore 
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where I t = hi + I 22 + I33 = n n M (a\ + a\ + a§)/5. 

Linearizing equations (A4)-(A5) including GR reaction alone, we have 
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and A eq = fl eq — Qm, where we have assumed that the GR reaction is a small correction. Linearizing 
equations for di, we obtain, for the bar-mode 
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With equation (A23), we then have 
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Again, this is only the frequency of oscillation in a frame with the complex angular frequency fl eq - 
The existence of an imaginary part in fl eq implies additional mode damping or growing. Similar to 
(A14), in a general frame with (real) angular velocity fi, the oscillation frequency is 

a ± (fl) = ±a - i " 25 Ql [5(Cl 4 M + 2tf M a 2 ) + af[ + 2 (fi - Q eg ) 

/_ n M \ 2K n Maj KtOm) 5 {A26) 
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Clearly, the Dcdekind-mode become secularly unstable when equation (A19) is satisfied, i.e., the 
gravitational radiation also induces an instability at the bifurcation point, but for a different mode as 
for viscosity. The growth time tqw is given by 



2n n Ma\ (fl M - cr a ) 
t gw = 25 ' > 



A. 4 Combined Effects of Viscosity and Gravitational Radiation 



The secular bar-mode instability tends to weaken or diminish when both viscosity and GR reaction 
are present. This is because viscous dissipation and GR reaction cause different modes to become 
unstable, as seen from §A.2 and §A.3. The combined effects of viscosity and GR reaction on the 
secular instability of incompressible Maclaurin spheroid were first noted by Lindblom & Detweiler 
(1977). Combining the results of §A.2 and §A.3, we find that the (complex) frequency for the bar 
mode in a rotating frame with angular velocity fMs given by: 
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For convenience, let us define dimensional viscous time constant and gravitational radiation time 
constant via 
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where the quantities in [• • •] are dimensionless and depend only on the eccentricity e (for compressible 
systems, there is also a dependence on n through R a /R), and 
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When t± < 0, the mode involved is stable, and |r±| then gives the corresponding damping time. In 
the limit of n = 0, the above results agree with that of Lindblom and Detweiler (1977). 

It is clear from the above expressions that the bar mode of a a compressible Maclaurin spheroid 
tends to be stablized by the combined effects of viscosity and GR reaction. To locate the region of 
stability, define 

Q^{l-elf'^- d \ (A33) 

"vis 
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where = 0.95289 and Q d = q n tt d /(irp) = 0.44022 refer to the dynamical stability limit. The 
Dedckind-mode is stable when 

QQ 4 d (}-^) ' (|~) (^m + (To) > (Am - to) 5 , {AM) 
and the Jacobi-mode is stable when 

Q{l d (}— J) ' (Jf) {&m - &o) < (fiM + Zof. (A35) 

Figure 9 depicts the critical T/|W| (where the secular instability first sets in) as a function of Q for 
n = 0, 0.5, 1, 1.5. For n = and Q = 1, the Maclaurin spheroid is stable all the way to the dynamical 
limit. 
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Figure Captions 



FIG. 1. — The secular instability growth time due to gravitational radiation as a function of the ratio 
[3 = T/\W\ for a neutron star modeled as a compressible Maclaurin spheroid, with M — 1.4M Q and 
R = 10 km. The solid line corresponds to the bar-mode (to — 2) instability for polytropic index 
n = 0, the dotted line for n = 0.5, and the short-dashed line for n = 1. The long-dashed line is for 
the m = 3 mode and n = 0. 

FIG. 2. Secular evolution tracks of a Riemann-S ellipsoid with n = 1 driven by gravitational 
radiation reaction. The energy of an ellipsoid as a function of the axis ratio a^lax is shown along 
various equilibrium sequences. Dcdekind-like sequences > 2|f2|) are shown on the left panel, 
Jacobi-like < 2|f2|) on the right panel. The thick solid curve corresponds to the Dcdckind and 
Jacobi sequences, the thick vertical line corresponds to the Maclaurin sequence, while the other lines 
correspond to constant-C sequences: C — C/(M 3 R ) 1 ^ 2 = —0.4 (dotted line), C — —0.32 (dashed 
line), and C — —0.25 (long dashed line). The solid round dot marks the point of bifurcation. The 
Dcdekind-like panel is divided into region I (fi > 0, C > 0) and region IV (fi < 0, £ > 0) by the 
Dedekind line (the solid curve), while the Jacobi-like panel is divided by the Jacobi line into region II 
(n > 0, C > 0) and region III (fi > 0, C < 0). 

FIG. 3. — The total energy radiated in gravity wave during the evolution of a star from a Maclaurin 
spheroid to a Dedekind ellipsoid as a function of the initial T/|W|. The solid line is for n = 0, the 
dotted line for n = 0.5, and the dashed line for n = 1. 

FIG. 4. — The gravitational waves emitted by a secularly unstable neutron star, evolving from a 
Maclaurin spheroid toward a Dedekind ellipsoid, (a) shows the gravitational wave amplitude as a 
function of the frequency; (b) the number of cycles spent near /, and (c) the timescale of the evolution. 
Here M = 1.4M Q , R = 10 km and D is the distance to the star. The solid lines correspond to the 
initial (3 = 0.24 and the dotted lines (3 — 0.2, both for n = 1; the dashed lines correspond to (3 = 0.24 
and the long-dashed lines (3 — 0.2, both for n = 0.5. 

FIG. 5. — (a) The maximum frequency f max of the gravitational wave radiated during the evolution 
of a neutron star from unstable Maclaurin spheroid to stable Dedekind ellipsoid as a function of the 
initial T/|W|. The solid line is for n = 0, the dotted line for n = 0.5, and the dashed line for n = 1. 
(b) Twice of the rotational frequency of the initial Maclaurin spheroid. Here M = 1.4M Q and R = 10 
km. 

FIG. 6. — Comparison between the characteristic amplitude h c of gravitational waves emitted during 
the secular evolution of a nonaxisymmetric neutron star and the rms noise h rms in LIGO's advanced 
detector. The heavy lines in the center of the figure correspond to the evolution shown in Fig. 4, while 
those in the right upper corner correspond to the evolution shown in Fig. 7 (but here only the solid 
and dotted lines are shown for clarity). 

FIG. 7. — Gravitational radiation from the evolution of Jacobi-like ellipsoids (regions II and III 
in Fig. 2). The quantities shown are labeled as in Fig. 4. The polytropic index is n — 1, while 
M — IAMq and R = 10 km. The solid lines and the dotted lines start from a Jacobi ellipsoid, end 
in a Maclaurin spheroid with (3 = 0.12 and (3 — 0.135. The short-dashed lines correspond to evolution 
toward a Maclaurin with (3 = 0.24, and the long-dashed lines to (3 = 0.02. 

FIG. 8. — Terminal frequency of the Jacobi-like evolution as a function of T/IT4 7 ! of the final Maclaurin 
spheroid. The vertical line marks the bifurcation point. 

FIG. 9. — The evolution of Dedekind-like ellipsoids (corresponding to region IV in Fig. 2) emitting 
gravitational waves which carry negative angular momentum, (a) shows the equilibrium energy, (b) 
the angular momentum, (c) the axis ratio, and (d) the gravitational wave amplitude, ft is the angular 
velocity of the ellipsoidal figure, so that the frequency of the wave is |fi|/7r. The solid lines correspond 
to C = C/(M 3 R ) 1 / 2 = —0.237, the terminal configuration is a Maclaurin spheroid with (3 = 0.1; the 
dotted lines correspond to C = —0.321, for which the final state is a Dedekind ellipsoid with £1 = 0. 
The polytropic index is n = 1. 

FIG. 10. — Competition of Viscosity and gravitational radiation on the stability of a Maclaurin 
spheroid for n = (solid line), n — 0.5 (dotted line), n = 1 (short-dashed line) and n = 1.5 (long- 
dashed line). 
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